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ABSTRACT 

A new vehicular traffic flow model based on a stochastic jump process in 
vehicle acceleration and braking is introduced. It is based on a master equa- 
tion for the single car probability density in space, velocity and acceleration 
with an additional vehicular chaos assumption and is derived via a Markovian 
ansatz for car pairs. This equation is analyzed using simple driver interaction 
models in the spatial homogeneous case. Velocity distributions in stochastic 
equilibrium, together with the car density dependence of their moments, i.e. 
mean velocity and scattering and the fundamental diagram are presented. 

1 INTRODUCTION 

Since the beginning of measurement of the vehicular flow on roads, it is well 
known that the dependencies of flow quantities show characteristic features 
mainly independent on individual drivers. Therefore a wide spread of the- 
oretical models are developed in the last decades. Today these models are 
roughly sorted into three categories: microscopic, macroscopic and meso- 
scopic models, see, e.g., [1, 2, 3, 4]. In microscopic modeling the acceleration 
of each car is given as a function of a subclass of all other considered kine- 



matic variables, e.g., [5, 6, 7, 8]. There, driver behavior is introduced using 
deterministic parameters or stochastic functions. Those models are often 
solved by computer simulations, e.g., [9, 10, 11, 12, 13]. The aim of macro- 
scopic modeling is building time propagation equations for the measured flow 
variables, i.e. car density, traffic density and mean velocity at each place of 
the road, e.g., [14, 15, 16, 17, 18]. These propagation equations are often 
of differential type using densities as continuous limits of the discrete traffic 
flow [19]. Vehicle number conservation on a road segment gives rise to the 
kinematic equation, which propagates the car density in analogy to the dif- 
ferential mass conservation law in fluid dynamics. The flow or mean velocity 
is modeled either by using a special ansatz or by constructing a second, dy- 
namic equation, like the Navier-Stokes ansatz in fluid dynamics. Up to now 
there is an open discussion of the applicability of dynamic equations to the 
description of traffic flow [20, 21, 22]. 

Another class of models follows the assumption that traffic flow is a 
stochastic process [23]. Here queuing theory methods or state probability 
distribution master equations are introduced, e.g., [24, 25, 26, 27]. Because 
nearly all models use information of the single vehicle behavior as input 
and produce results of the whole traffic flow, this model type is often called 
mesoscopic. The mesoscopic ansatz was applied to single vehicle states, e.g., 
[28, 29, 30, 18, 31, 32], as well as for vehicle clusters [33, 34]. Their deriva- 
tion is analog to the equations in gas kinetic transport theory, especially the 
Boltzmann transport equation, though there is no momentum and energy 
conservation in the interaction between two vehicles. Models based on single 
car state distributions are widely called kinetic models in traffic flow. 

There are two principle binary interaction types in kinetic models. The 
acceleration of a vehicle due to the lack of leading vehicles and the braking 
due to a leading vehicle. Overtaking normally is introduced via a special 
passing probabflity depending on the car density, [4], or using multi lane 
coupling [35]. Though acceleration often is modeled as a mean- value func- 
tion [29, 27], the de acceleration interaction is a jump process in the velocity 
variable, i.e. vehicles have infinite braking strength with vanishing durance. 
At the end of a de acceleration interaction the following car has the same 
velocity as the leading one. Kinematic variables of the leading car remain 
unchanged. Due to the velocity jump process kinetic models include the im- 
plicit assumption that the interaction durance is much smaller than the mean 
free driving time between two interactions. Additionally a 'vehicular chaos 
assumption' must be fulfilled [29, 31] . Car following experiments on road seg- 
ments show that velocity changes due to interactions with leading cars take 
several seconds typically [36], where moderate acceleration/braking values in 
low or medium dense traffic are assumed. Under these conditions during an 
interaction another interaction can occur, producing multiple interactions. 
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which are not included in those kinetic models. This effect increases with in- 
creasing car density because the overtaking probabihty decreases. Therefore 
kinetic models are often applied to partly constrained traffic [2] . 

There are several possibilities to enhance kinetic models into the higher 
car density regime. One of them is to introduce a continuous velocity change 
due to a mean acceleration and a mean velocity scattering. Then the single 
car state probabihty density is described by a Vlasov-Fokker-Planck type 
equation [27, 37] . Another possibility is to introduce the acceleration variable 
into the underlying stochastic process to be the main driver control variable 
of the interaction, changed in a discontinuous manner. The second ansatz is 
followed in this paper. 

The traffic flow model introduced in the following sections is a master 
equation approach for the single car state probability density function simi- 
lar to kinetic models. The main difference to these models is the jump process 
in the acceleration variable, whereas the velocity changes continuously and 
depends on the acceleration in a deterministic way. As is shown in car fol- 
lowing measurements acceleration changes take in the order of one second 
[36] . So except in extreme cases the assumption of binary interactions seems 
to be valid even at higher densities. Acceleration changes depend not only 
on vehicle kinematics but also on the driver skill and behavior [38], which 
introduces an independent component to the model and therefore justify the 
vehicular chaos assumption. 

The model is based on a car following ansatz, using kinematic variables of 
vehicle pairs, building a stochastic vector, propagating in time as a Markov 
process following the Feller-Kolmogorov equation [39]. This equation is re- 
duced into an equation for calculating the single car state probability den- 
sity by generalization of the vehicular chaos theorem used in kinetic models 
[31, 32]. For model evaluation a simple interaction function between the cars 
of a leading car pair is constructed and applied to a homogeneous traffic flow 
without overtaking. The influence of different interaction functions to the 
velocity distribution and the existence of stochastic equilibria together with 
the car density dependence of the typical flow features are discussed. 

Note that a lot of material used in the following section is standard in 
mathematics of stochastic processes, although some repeating is important 
to clarify the underlying model assumptions. 

2 DERIVATION OF THE MODEL EQUATION 

This section is divided into two subsections. First, the stochastic process is 
defined and second, using a vehicular chaos assumption, a master equation 
for the state probability density of single cars is constructed. 
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2.1 Definition of tlie Stocliastic Process 



For a given car pair at time t the leading car is located at Xt with velocity Vt 
and acceleration at and the following car is located at Xt < Xt with velocity Vt 
and acceleration at defining a state vector = {xt, Vt, at, Xt, Vt, at). As time 
increases from t to t + r, where r is infinitesimal small, Xf, Vt, Xf, Vf changes 
continously to Xt+r, Vt+r, Xt+r, Vt+r due to the kinematic laws. Restricting 
the model to binary interactions, where the leading car is not influenced 
by the following one, i.e. at+r — o^t, the acceleration of the following car 
at+r changes discontinuously depending on the state vector y^. The jump 
hight then is given by e,- = at+r — dt- Because in general is distributed 
stochastically with probability density Pj(er|yt) the state vector is a random 
vector Yt with representation y^ at time t, therefore following a Markov 
process. The transfer probability density p from state = y^ to state Yt+r 
is given by 

p(Yt+r|yt) = 

5{xt+r - Xt- VtT - -atT'^)5{vt+r - Vt - atT)5{xt+r - Xt - VtT - i^atr'^) ■ 
S{vt+T -Vt- atT)5{at+r - at) ■ {5{at+r - at) + Pj{er ^ 0|yt)) (1) 

where b{z) is the Dirac function, [40], representing the deterministic, contin- 
uous parts of the change process. During r there can be an acceleration jump 
due to density or no change, represented by an 5-function. Because both 
events arc incongruous, they appear in Eq. 1 as a sum in brackets. Note that 
Pj can be written conditioned, Pjie^ ^ 0|yj) = pj[er\er 7^ 0, yt)-Pj{^r 7^ 0|yt), 
where pj is the density of a change of height under the condition that a 
jump occurs and Pj is the probability of a jump. 

Defining the state probability density of a leading car pair f^iyt) — 
/2(y,t) at time t, the time propagation of this quantity is described by the 
Feller-Kolmogorov equation [39] 

^/2(y, t)= f (g(z ^ y, t)/2(z, t) - q{y ^ z, t)h{y, t)) (2) 

with initial condition /2(x, i = 0). Here q is the density of the transfer rate 
between states and is connected to the transfer probability density p, Eq.l, 
as 

dp 



(3) 

=0+ 

Normally the components of the state vector z are defined in some given 
intervals producing restrictions to the integration area Q. Here Q is expanded 
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to the whole state vector space by continuing the state probabihty density 
/2 zero valued outside of the restricted area. So possible source terms are 
shifted into possible boundary conditions. The influence of the restrictions 
will be discussed below. 

Now Eqs. 1 and 3 are inserted into Eq. 2 and some algebra with S- 
functions and their derivatives are done [40]. Together with the total transfer 
rate Q, defined by 



Q(y,t) = q{y^ z,t)dz 

Jzen dr 

and the calculated integral 



(4) 

T=0+ 



. ,w / df2 df2 _df2 _df2 

q{z^ y,t)f2{z,t)dz^ -V- a- v— - a— 

xeu ox ov ox Ov 



L 

+ / Q{x, V, a', X, V, a, t) ■ Pj{eo = a — a'\eo ^ 0, x, v, a', x, v, a, t) ■ 

Ja'eR 

f2{x, V, a' , X, V, d, t)da' , (5) 

where z = {x', v', a', x', v', d') and eo = £7^=0+ is used, the Feller-Kolmogorov 
equation reads as 

«9/2 , 0/2 df2 _df2 _df2 
ot Ox Ov Ox Ov 

= / Q{x,v,a' ,x,v,d,t) ■ pj{€o — a — a'\€o 0,x,v,a',x,v,d,t) ■ 

f2{x, V, a', X, V, d, t)da! . (6) 

The function arguments on the left side are suppressed for simplification. 
There are two important features not included in Eq. 6 until yet. First, the 
state of the leading car (x, d) only changes continously. So it has not the 
same character than a following car, where acceleration jumps can occur. 
Second, the equation does not include the property x < x for a leading car 
pair. Both features will be introduced in the next section via a vehicular 
chaos theorem, splitting /2 partially in a product of single car probability 
densities. 



2.2 Vehicular Chaos Assumption and the Master Equation 

Vehicular chaos assumptions are often used in Boltzmann type mesoscopic 
traffic flow models. The main idea is to decouple /2 into a product of two 
single car probability densities neglecting correlations. In this paper the ideas 
of R. Wegener and A. Klar [32] and P. Nelson [31] are essentially used. 
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The car pair state density /2 can be written as a product of conditional 
densities as 



f2{x, V, a, X, V, a, t) — li{x, v, a, t) ■ l2{x\x, v, a, t) ■ ^3(1', a\x, x, v, a, t) , (7) 

where the first factor li is identified with the single car probability density 
f{x,v,a,t) and the second factor I2 is the conditioned distance density D 
introducing the distance variable h — x — x instead of i.e. 

l2{x\x^ V, a, t)dx = D{h\x, v, a, t)dh . (8) 

The third factor describes the behavior of the leading car depending on 
the state of the following one and therefore should be a function of /. The 
vehicular chaos ansatz assumes that the velocity v and the acceleration d of 
the leading car is independent on the state of its following one resulting in a 
special form of Z3, 

h{v, d\x, X, V, a, t) = k{v, d\x, t) = "^^^^^'^^ ' (9) 

with definition 

T{x,t) = I f{x,v,d,t)dvdd. (10) 

J v,a 

This special form of /2 ensures both features not included in Eq. 6. Each 
leading car now is distributed as v, a, t) and therefore is also a following 
one in another leading car pair. The distance density D is zero valued for 
distances h < /imin, which is equal to the mean car length. It is a behav- 
ioral function and therefore does not depend explicit on a special place x 
at time t. However, there is an additional dependence of the distance on 
the overall flow, mainly represented by its mean values like mean velocity, 
scattering of velocity and acceleration, etc.. Typically these are lower mo- 
ments of the single car state probability / and so space and time dependent. 
They arc included into D as elements of a moment vector mf(x,t) of /, i.e. 
D = D{h\v,a.mf{x.t)). Note that this moment vector can also include fur- 
ther flow dependent functions like the mean distance, which are not direct 
moments of / ( see Eq. 13 below). 

Inserting this vehicular chaos assumption in Eq. 6, integrating over a, v 
and X or h and bearing in mind that the single car state density / vanishes 
at the infinite boundaries of v and x continuing / like /2, the resulting time 
propagation equation for / is given by 

df df df f f{x + h,v,d,t) 



— — \-v— — \- a— = I dhdv dada ■ {1 — Po) ■ 

Ot ox OV Jh.v.d.a' 



J^{x + h,t) 
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{a{a\x, h, V, v, a' , a, t)Q{x, h, v, v, a', a, t)D{h\v, a', inf{x, t))f{x, v, a' , t) 
— a{a'\x, h, V, v, a, a, t)Q{x, h, v, v, a, a, t)D{h\v, a, mf{x, t))f{x, v, a, t)} 



introducing 

Pj{eo — a — a' — \eo ^ 0, x, v, a', x, v, a, t) = (1 — Po)cr{a\x, v, a', h, v, a, t) . 



Here a simple ansatz for overtaking, often assumed in kinetic traffic models [4] 
is used. It is based on the overtaking probability density Pq, which depends 
on the kinematic variables of the car pair and other vehicles nearby. Because 
of its mathematical complexity, in literature this function often is simplified 
on the local car density only. For a more precise modeling of overtaking 
effects a multilane model has to be constructed in the same manner as is 
done elsewhere [35]. 

The car density K at place x and time t is given by the inverse of the 
mean distance H at the same place and time, i.e. K{x,t) = 1/H{x,t). This 
unconditioned mean value can be calculated using the vehicular chaos ansatz 
resulting in 



Together with Eq. 13, Eq. 11 describes the whole traffic dynamic. Note that 
if D is independent on the state of the car (v, a), Eq. 13 is always fulfilled. 

Eq. 11 has the same structure as those, used in Enskog theory of dense 
gases [41] or other kinetic traffic flow models in literature [4, 27] with a spa- 
tial correlation via D. It couples microscopic car pair interactions deflned 
by a, Q and D with macroscopic traffic flow quantities, which are moments 
of the state probabihty density /. The main differences to other traffic flow 
models based on kinetic type equations are the additional, independent ac- 
celeration variable, free functions describing the car pair interaction and the 
model derivation from a Markov process. In contrast to the often used mean 
acceleration ansatz, here the acceleration is used as the control variable of 
the driver process as in reahty [29]. 

Note that using the same derivation technique a model based on a con- 
tinuous stochastic change of acceleration can be derived, which is similar to 
the Fokker-Planck-equation framework [42], but in the acceleration variable. 
There, in contrast to the jump process the time sequence of acceleration 
changes due to interactions has to be known as input, making the model 



(11) 



(12) 





f{x,v,a,t) 
J^{x,t) 



dh dv da . (13) 
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more complex and less applicable [43]. 

For solving Eq. 11 an additional initial condition and spatial and velocity 
boundary conditions must be specified. The velocities of all cars have to be 
non negative, i.e. v > 0, and the acceleration a is bounded Omin < a < Omax, 
where amin < is the maximal braking value. If a speed limit is defined 
on a given road segment, w{x), then all velocities are lower than this limit 

V < w{x). Normally the spatial coordinate x lies between a minimal and 
a maximal value, where in flow and out flow conditions must be specifled. 
These restrictions also influences the construction of the interaction functions 
prohibiting values outside of the allowed intervals. Concrete interaction ex- 
amples for model testing will be given below. 

Driver measurements show that interaction decisions are mainly done by 
the own driving state v, a and the distance h and roughly the relative velocity 

V — V between the cars in a leading car pair [44, 45, 3]; the acceleration of 
the leading car is not taken into account. 

Assuming a spatial homogeneous road segment and the arguments dis- 
cussed above, Eq. 11 is reduced to a Boltzmann type equation 

df df f _ 

— + a— = / dv da' ■ f{v, t) ■ 

at OV Jv,a' 

{T,{a\v, V, a , mf{t))f{v, a , t) — T,{a\v, v, a, m/(t))/(t', a, t)} , (14) 
with f{v, t) — f{v, a, t) da and 

'E{a\v,v,a' ,nif{t)) = / a{a\h,v,v,a')Q{h,v,v,a')D{h\v,a' ,m.f{t)) dh 

*' ^min 

(15) 

is the weighted interaction density. Before the complex spatial dynamics of 
a traffic flow is studied in a model, the typical phase diagrams in stochastic 
equilibrium are calculated and proofed to be in agreement with measured 
data. In contrast to gas kinetic, where the universal Maxwell-Boltzmann 
distribution exits, in traffic flow it depends on chosen interaction functions, 
because there are no Stofiinvarianten in the process. The simplest way to 
calculate solutions in stochastic equilibrium is starting from an homogeneous 
initial condition and following its time propagation assuming a homogeneous 
road segment using Eq. 14. Naturally the real process can not assumed to 
stay homogeneous during time propagation, because fluctuations in spatial 
dynamics occur, vanishing when reaching the equilibrium. 
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3 SOME REMARKS ON THE MACROSCOPIC MODEL BEHAVIOR 



In contrast to stochastic equilibrium, the spatial dynamics in traffic flow 
models are often studied using macroscopic flow equations calculated from 
the base equation using the moment method. Although this paper deals with 
direct solutions in the equilibrium, some preliminary remarks on the macro- 
scopic model behavior are given to facilitate a comparison to other models. 
Additionally for a better understanding of the model structure and compar- 
ison to other kinetic models the equation for the reduced state probability 
density f{x, v, t) — f{x, v, a, t) da is discussed. This equation is calculated 
integrating Eq. 11 with reference to a. At flrst sight, the result 

df df dAf , , 

seems to be a Vlasov-like equation with a mean acceleration fleld A. This 
function is calculated by the flrst moment of Eq. 11, 

dA] dAf dSf f 

where S = J^a^JXxjV, a, t) da and the operator J' is the right side of Eq. 11. 
So A depends not only on the kinematic state of the car, but also on the 
higher moments of / especially on the scattering S. Now S can be calculated 
by the next moment of /, etc.. For practical use a closure condition must 
be introduced resulting in a system of differential equations. As mentioned 
above in literature the Vlasov-Fokkcr-Planck ansatz is often used to model 
acceleration behavior in a mean value sense with velocity scattering. In 
contrast to those models, here an explicit velocity diffusion term can not be 
generated straight forward, but it is implicit included into the function A via 
Eq. 17. Up to now it is an open question, under which conditions an explicit 
diffusion term can be constructed in a consistent way. 

Macroscopic equations are produced by velocity and acceleration mo- 
ments only depending on x and t. Integrating Eq. 16 with reference to v, 
using Eq. 10 yields the continuity equation 

The equation for the mean velocity V{x, t) results analogous calculating the 
flrst velocity moment 

dV ^dV _ ^_de _e dJ" 
dt dx dx T dx^ 



9 



with mean acceleration B{x, t) = A{v, x, t)dv and velocity variance 0(a;, t) = 
Iv ^^/(^) t)dv — V'^{x, t). B and © can be calculated by the next moments 
of Eqs. 16 and 17. Mean acceleration B and mean velocity V are coupled by 
the mean covariance and the acceleration scattering in the second velocity 
moment of Eq. 17, etc.. The structure of both first order Eqs. 18, 19 are 
similar to those known from miscellaneous literature models [46]. There B is 
specified by an ansatz in contrast to this model. Up to now closure conditions 
needed for practical use of the macroscopic equations are not developed and 
therefore have to be done in a future work. 

4 EXPLICIT SOLUTIONS IN STOCHASTIC EQUILIBRIUM 

In this section the model equation is solved for special interaction functions. 
The aim is to show that under justifiable input the result also seems to be 
applicable. Additionally the section shows the existence of equilibria for 
these exemplary interactions. First, a simple relative velocity dependent in- 
teraction, which can be found in low density traffic fiow is developed and its 
solution behavior is discussed. Second, a simple distance threshold interac- 
tion is introduced and the car density dependence of the flow quantities is 
shown. 

4.1 Relative Velocity Dependent Interactions in a Low Density Traffic Flow 

At low car densities, the occurrence of an interaction together with the inter- 
action strength are not depending on the distance between the cars, but on 
the change of the distance, i.e. the relative velocity mainly [45]. The driver of 
the following car has nearly no quantitative estimates on the kinematic state 
of the leading car and therefore decides his new acceleration value by his own 
behavior. Only the sign of the new acceleration value, braking or accelerat- 
ing, depends on the sign of the relative velocity. The simplest approximation 
of this behavior is a two value acceleration change, given by 

a{a\v — v) = Q{v — v)6{a — ai) + Q{v — v)6{a — 02) , (20) 

where Q{z) is the Heaviside step function and ai > is the fixed new ac- 
celeration value and 02 < is the fixed de acceleration value. Assuming 
this relative velocity dependent interaction, the interaction rate Q is only a 
function oi v — v. In this model all driver behave similar and therefore have 
nearly the same velocities. So Q can be Taylor expanded for small relative 
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velocity values up to first order as 

Q{v-v)^Qo + ro\v-v\ . (21) 

In the following, two extreme cases will be considered: the relative velocity 
rate Qo — 0, which in gas kinetic is called hard sphere rate, and the constant 
rate Tq — 0, which in gas kinetic is called a Maxwell rate. Because in this 
section a homogeneous traffic flow based on Eq. 14 is supposed, Eq. 15 is 
given by 

S(a|-u — v) = Q{v — v) • a{a\v — v) , (22) 

which is independent on the chosen concrete distance density D. Therefore 
the velocity distribution is also independent on distance features. To produce 
car density dependencies, the solution f{v, a, t) must be inserted into Eq. 13 
and a concrete D must be defined there. 

Note that this interaction profile does not ensure positive velocities. So it 
can only be applied, if the mean velocity of the flow is large against its veloc- 
ity variance, which lets the probability for self-defeating negative velocities 
vanish approximately. In the next subsection it is shown, how such negative 
velocities can be avoided using velocity boundary conditions. 

Assuming an initial condition, where all cars have one of the two possible 
acceleration values, the state probability density can be written as 

f{v, a, t) = /i {v, t)S{a - ai) + Mv, t)S{a - a^) . (23) 

/i and /2 are the state densities of the accelerating and the braking cars 
resp.. Inserting Eq. 22 and Eq. 23 into Eq. 14 results in two equations 

f2(v,t) [_ Q(v-v)f(v,t)dv 
Mv,t) [_ Q{v-v)f{v,t)dv, 

Jv<v 

h{v,t) [_ Q{v-v)f{v,t)dv 
f2{v,t) [_ Q{v-v)f{v,t)dv. (24) 

Jv>v 

From both equations it can be shown that the solutions after a sufficient long 
period of time t are given asymptotically by 

Mv,t)r^^f,^{v-A-t) , i = l,2, (25) 

where A = {ai + a2)/2 is the mean acceleration and /cq is the equilibrium 
solution of the same system, Eq. 24, whereas Oi is replaced by (ai — a2)/2 



dt dv 



dt dv 
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and a2 is replaced by (02 — ai)/2, . Note that for these new acceleration 
values the mean value A vanishes and so a stochastic equilibrium exists. 

Starting from an initial condition for /j, after some time the shape of the 
functions are changed into the shape of fcq, which then propagates along the 
velocity axis with mean acceleration (ai+a2)/2. Figure 1 shows an example of 
this time evolution starting with two Gaussian shape densities and resulting 
in a moving one given by Eq. 31 below. It is calculated numerically by 
discretizing Eq. 24 using a standard exphcit discretizing scheme for the left 
hand side and the standard trapezoidal rule for the integrals on the right 
hand side [43]. 

To calculate the equilibrium solution, first ai = —02 = Oq is set. For 
dfi/dt — and using feq{v) = fi{v) + f2{v), which is calculated directly from 
Eq. 23, the sum of both Eqs. 24 results in an equation for the equilibrium 
velocity density 

df^ ^ 4M ( r _ v)U^)dv - r Q{v - v)Uv)dv) (26) 

with additional condition f(.{v)dv = 1. The difference of both equations 
results in fi{v) = f2{v) and so fi{v) = feq{v)/2 for i = 1,2. The general 
solution, Eq. 23, is 

N ? / J{a-ao) + S{a + ao) 
f{v, a) = /eq(v) . (27) 

For rate Q — rQ\v — v\ Eq. 26 is reduced to 

%M^Mi(^_,) (28) 
dv Oq 

with a Gaussian probability density as solution 

/e» = -^-ie-(^-W. (29) 
V27r CTy 



Here the mean velocity y is a free parameter and ay — ^Jao/r\,. In the case 
of a constant rate Q — Qo — l/T, Eq. 26 is reduced to the system 

Taking the derivative of the second equation, inserting the first and then the 
second again, results in s"{v) = [s(f )^]'/2Tao, which can be solved using the 
method of variable separation. After inserting s{v) in the second Eq. 30, the 
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solution is 

with free parameter V again and — aoTn/y/S. Figure 2 shows both ve- 
locity distributions in standardized form. The deviations are small and the 
Maxwell type density is slightly less peaked. Both distributions agree quali- 
tatively to measured distributions and other theoretical work [47, 27]. Note 
that the equilibrium velocity distribution seems not to be strong influenced 
by the functional form of Q in this traffic flow regime. This result is the 
same as is found for thermalization processes governed by the homogeneous 
Boltzmann equation in kinetic gas theory [48]. There also the shape of the 
energy distribution during an equilibration process is nearly independent on 
the interaction rate; the interaction rate influences mainly the timescale of 
thermalization . 



4.2 A Simple Distance Threshold Interaction Model 

In this subsection a distance oriented interaction model is introduced, to 
study the car density behavior of the equilibrium distribution. Here a simple 
distance interaction threshold H{v) = + well established in literature 
[49] is deflned. Only if the distance between the cars in a leading car pair is 
equal to H{v) an acceleration change is done by the following car. The driver 
of the following car checks this condition with a constant rate Q — Qq — 1/T. 
Like in the last subsection, here again a two value acceleration change profile 
is used for simplicity, i.e. 

a{a\h, v) = e{h - H{v))S{a - oq) + e{H{v) - h)S{a + oq) . (32) 

If the distance h is larger than H{v), then the following car accelerates with 
a = qq, otherwise it de accelerates with a = —gq. Because a depends on h, 
a concrete distance correlation D is needed. As discussed in hterature the 
exponential ansatz, unconditioned on {v, a), seems to be the right choice for 
moderate car densities [31] 

where H is the mean distance between the cars in a leading car pair and 
therefore is equal to /K. Note that this ansatz fulfills Eq. 13. Inserting 
this interaction profile into Eqs. 14, 15 analog to the last section, esp. using 
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Eq. 27, the equation for the equihbrium velocity distribution is 



df (v) f iv'\ r^iv) 
uj^Kj^l ^ ^^^(1 - 2P{h < H{v))) , with P{h < H{v)) = / D{h) dh . 
dv aoT Jhmir, 

(34) 

The interaction profile does not avoid the unreahstic negative velocities. For 

V < there is H(v) < hmin and therefore P{h < H(v)) — 0. The solution 
of Eq, 34 then is an exponential function decreasing with decreasing v. For 

V > Eq. 34 together with Eq. 33 can be solved analytically. Both solutions 

V < and v > are matched at v = 0, because there the function is 
continuous. After some ponderous but straight forward calculations including 
the unity norm, the probability density is given by 

^'''^''^ " aor(e-2/3 + /3(2/?)-/37(/3,2/?)) ' ^^^^ 

where 7(0;, y) is the incomplete gamma function and (3 = {H — hjam) / {cxaoT) 
the scaling parameter of the solution. It compares the mean distance between 
two cars with the distance change during T. From this solution analytic ex- 
pressions for the mean velocity and the velocity scattering can be obtained. 
They are not included here, because of their mathematical complexity. The 
intrinsic model error P{v < 0) = /^^^ /eq(^) dv vanishes exponentially pro- 
portional to {e/2)~f^/^/P for large P and so for small car densities. 



4.3 The Introduction of Velocity Boundary Conditions 



To eliminate negative velocity probabilities absorbing boundary conditions 
at v = in analogy to the Bosc gas state condensation can be constructed. 
If a braking car reaches f = 0, it changes independent on the state of the 
leading car his acceleration value to a = 0. Introducing a maximum free flow 
velocity w, the same ansatz is possible there. If a car reaches v — w, its 
acceleration value is set to zero. Using both boundary conditions the new 
state density f{v,a,t) is calculated from 

f{v, a, t) = f{v, a, t)x(0 < v < w) 

/O _ rco _ 

f{v, t) dv 5{v)6{a) + / f{v, t) dv 6{v - w)6{a) . (36) 
-00 Jw 

x(0 < v < w) is the index function, which equals one ior < v < w and 
zero elsewhere. This boundary condition influences the behavior of velocity 
distributions especially in the free flow regime and near v = 0. It also changes 
the equilibrium behavior together with the existence statements presented in 
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the last subsection. So Eq. 24, its solution 25 together with the boundary 
condition results now in an equilibrium solution feq{v, a) = 6{a)6{v) for A < 
and feqiv, a) — 5{a)5{v — w) for ^4 > 0. 

Applying the boundary condition to the distance threshold model, sub- 
section 4.2, /cq(f,a) is given by Eq. 27 together with Eq. 35. Integrating 
Eq. 36 with reference to a results the new velocity distribution fcq{v), which 
is shown in figure 3 for different car density values K = 1/H and two dif- 
ferent acceleration values oq. The distributions show the typical shape as 
mentioned above and their situation depends strongly on the car density. 
Note that Qq, which is also equal to the acceleration scattering (acceleration 
noise, ACN) in this simple interaction model mainly influences the width or 
scattering of the distributions. 

The car density dependence of the mean velocity V — Jq vfeq{v)dv and 
the velocity scattering calculated from — v'^feq{v)dv — V'^ are shown 
in the figures 4 and 5. At low car densities the mean velocity equals the free 
flow velocity w due to the boundary condition. With increasing car density 
the mean velocity shows the typical decrease. Because for a = 1.8s an inter- 
action occurs at larger threshold distances than for a — 1.2s, the decrease 
occurs at lower car densities. Velocity scattering results are much more sel- 
dom in literature than mean value results. The model shows a steep increase 
at moderate densities, where the transition from free flow to interaction ori- 
ented flow occurs. Then with increasing car density the velocity scattering 
diminishes. Other threshold oriented kinetic models show the similar behav- 
ior [49] . The peak in the scattering shifts to higher car density values when 
a decreases, which corresponds to the begin of the mean velocity decrease in 
flg. 4. The fundamental diagram, flgure 6 shows its typical behavior. The 
interaction oriented part of the curve depends on the a value. Let the lower 
a-value is identified to a more aggressive driver and the higher a-value to 
a more conservative one, a traffic fiow of a mixture of both produce traffic 
densities in the area between both curves. This could be a first guess to mul- 
tivalued traffic flows in the framework of this model as discussed in literature 
[50, 37]. 

In gas kinetic models interactions produce velocity jumps based on inflnite 
acceleration and de acceleration change values. In the model discussed in the 

last two subsections this corresponds to the case of Oq — > cxo. Applying this 
limit operation to the boundary condition for all three solutions Eq. 29, 31 
and 35, together with Eq. 27, a limit probability density independent on the 
concrete interaction proflle but based on velocity jumps can be calculated 

Uv, a) = l^^5{v) + ^S{v - w)^ ■ 5ia) . (37) 
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The velocity part of this density is well known in kinetic traffic flow theory 
and used there as an argument for the existence of bimodal velocity distri- 
butions [31, 37]. There, [31], the factor 1/2 is replaced by a more complex 
function of the distance probability at v = w, but this difference is mainly 
due differences in the interaction model. The bimodal limit distribution in 
principal can be derived from a more universal qualitative argument. As long 
as the velocity scattering depends on the acceleration scattering in a mono- 
tone increasing way, the limit case together with the boundary conditions 
always produce this kind of solution with perhaps different prefactors. Its 
practical applicability seems to be problematic, because typical peak accel- 
eration scattering values are in the order of Im/s^, far away from this limit 
case [51]. 

5 CONCLUSIONS 

(1) In this paper a vehicular traffic flow model based on a stochastic car 
pair process with acceleration jumps is introduced. The model equa- 
tion for the single car distribution function in the spatial, velocity and 
acceleration variable is of Enskog type using a vehicular chaos assump- 
tion. The model input is based on car following information, especially 
for the interaction rate, the interaction strength and the distance distri- 
bution. The behavior of solutions for different interaction models and 
macroscopic moment equations are discussed. Velocity distributions in 
stochastic equilibrium arc calculated for special interaction cases. In 
contrast to Boltzmann theory, the existence and shape of an equilibrium 
solution of the model depends on the construction of the interaction 
functions. 

(2) The spatial homogeneous equation for a discrete two value acceleration 
change, depending only on the relative velocity, is discussed. There an 
equilibrium solution exists exclusive, if the mean acceleration vanishes. 
For this case velocity distributions for constant rate and relative veloc- 
ity dependent rate arc calculated analytically, showing qualitative the 
same behavior as is found in measurements or other traffic flow mod- 
els. With increasing acceleration scattering the velocity distribution 
broadens, i.e. the velocity scattering increases also. 

(3) A simple distance threshold oriented interaction model with constant 
interaction rate is discussed in stochastic equilibrium. To prohibit un- 
realistic negative velocities and to introduce a maximum free flow ve- 
locity, absorbing boundary conditions in analogy to the Bose gas state 



16 



condensation are introduced. The car density dependence of the veloc- 
ity distributions, the mean velocity and velocity scattering are shown 
for special interaction values. The fundamental diagram agrees quali- 
tatively in shape with other models. 

(4) In a next step the interaction profile must be enhanced to be more 
realistic. Especially the acceleration noise, here given by a constant 
value, together with the car density dependence of the acceleration 
distribution must be considered. Under these conditions the model can 
be solved only by computer simulation techniques. So a stable method 
to solve the master equation must be developed. Further work on this 
field is under way. 
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FIGURE CAPTIONS 



Fig. 1 : Time evolution of the partial distribution functions /i (broken lines) 
and /2 (solid lines) using Eqs. 24 for Maxwell interaction with param- 
eters T — 2s, oi = 0.4m/s^ and 02 = — 0.2m/s^. The initial conditions 
at t = are normal distributions with mean V — 28m/s and variance 
O.lmVs^. 

Fig. 2: Standardized velocity distributions in stochastic equilibrium based on 
Eq. 26. The brocken line describes the Maxwell type interaction rate, 
Eq. 31, where the solid line is the relative velocity type interaction rate, 
Eq. 29. 

Fig. 3: Equilibrium velocity distributions for different car density values K 
described in subsection 4.3 with oq = 0.05m/s^ (solid line) and oq = 
0.2m/s^ (broken line). The following fixed parameter values are used: 
T — 2.5s, a =1.8s, hmin — 6.5m. 

Fig. 4: The car density K dependence of the mean velocity V for two differ- 
ent threshold slopes as described in subsection 4.3. The following fixed 
parameter values are used: T — 2.5s, Oq =0.2m/s^, /imin = 6.5m. 

Fig. 5: The car density K dependence of the velocity scattering for two 
different threshold slopes as described in subsection 4.3. The following 
fixed parameter values are used: T — 2.5s, ao =0.2m/s^, hmin — 6.5m. 

Fig. 6: The car density K dependence of the traffic density K ■ V for two 
different threshold slopes as described in subsection 4.3. The following 
fixed parameter values are used: T — 2.5s, ao =0.2m/s^, h^am — 6.5m. 
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